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(Dated: February 1, 2008) 

Several proposed schemes for the physical realization of a quantum computer consist of qubits 
arranged in a cellular array. In the quantum circuit model of quantum computation, an often 
complex series of two-qubit gate operations is required between arbitrarily distant pairs of lattice 
qubits. An alternative model of quantum computation based on quantum cellular automata (QCA) 
requires only homogeneous local interactions that can be implemented in parallel. This would be a 
huge simplification in an actual experiment. We find some minimal physical requirements for the 
construction of unitary QCA in a 1 dimensional Ising spin chain and demonstrate optimal pulse 
sequences for information transport and entanglement distribution. We also introduce the theory of 
non-unitary QCA and show by example that non-unitary rules can generate environment assisted 
entanglement. 

PACS numbers: 03.67.-a, 03.67.Mn 



I. INTRODUCTION 

Much progress has been made recently in developing 
architectures that can support quantum information pro- 
cessing (QIP). The key result on the universality of quan- 
tum computers (QC) is that given the ability to imple- 
ment single qubit and two-qubit gates in a network of 
connected qubits any computation or quantum simula- 
tion can be implemented. However, some systems are 
better adapted to implement certain QIP tasks than oth- 
ers. In particular, lattice based systems with regularly 
arranged qubits interacting with nearest neighbors such 
as a neutral atom optical lattice, quantum dot arrays, 
and Phosphorus embedded Silicon offer several advan- 
tages in terms of reconfigurability and scalability. For 
example, optical lattices with nearest neighbor tunneling 
couplings have been shown to be a promising platform to 
simulate many body Hamiltonians [l[ . Generally, lattices 
are well suited to perform parallel computation protocols 
such as entanglement distribution |2| and entanglement 
swapping y[. 

The natural facility of these systems invites study of 
other models of computation that take advantage of the 
lattice architecture. Perhaps the most relevant computa- 
tional model in classical systems is a cellular automaton. 
The essential idea behind cellular automata (CA) is to 
make use of simple local rules uniformly applied across a 
lattice of cells to generate complex dynamics. Depending 
on the initial state of the system and the underlying rule, 
long range spatial and temporal correlations can develop 
resulting in complex behavior. Classical CA can simu- 
late a wide range of complex physical phenomena includ- 
ing fluid dynamics, nonlinear diffusion, percolation, and 
phase transitions in many body systems Formally, 
a CA is termed complex if it evolves in a manner that 
in some sense is computationally irreducible, meaning it 
cannot be predicted with a compactly written equation 
H . A number of C A rules have been shown to be compu- 



tationally universal, in the sense that they can emulate 
a universal Turing machine Q. 

The extension of the cellular automaton concept to 
quantum systems is fairly straightforward, though as we 
will show, requires a slight modification of the classical 
CA procedure for implementing the local rules. For a 
two-state CA, which will be the focus of this paper, each 
cell in a quantum cellular automata (QCA) corresponds 
to a qubit that can be in a superposition of states |0) and 
|1) and the local rule is carried out via a unitary gate op- 
eration on each neighborhood. The essential new feature 
in a QCA that makes it distinct from its classical coun- 
terpart is that nonlocal correlations can develop between 
cells resulting in the spread of entanglement throughout 
the system. This property of QCA will be of central im- 
portance in this paper. 

Our motivation for studying QCA is to explore the 
power of low computational depth circuits applied in uni- 
form across a system to produce complex quantum dy- 
namics. This is in marked contrast to the typical QC 
approach, where a complex sequence of logic gates act- 
ing on distributed qubits in the computer is carried out 
in a serial fashion in order to produce the desired output 
of a specific computation. Most previous work on QCA 
has focused on mapping such systems to the QC circuit 
model |6|,Lj|. Additionally, there have been investigations 
of quantum lattice gas automata (QLG) for simulations 
of the Dirac equation in ID 8] and for topological compu- 
tation 9] . Recently there was an experimental realization 
in liquid state NMR of a QLG algorithm to solve the ID 
diffusion equation 10]. We propose using ID QCA to ex- 
plore complex quantum correlations generated by simple 
rules applied over small neighborhoods. Characterizing 
multi-particle entanglement is a field of active research 
both for its potential use in QIP and in the study of non- 
locality in physics. QCA can offer a unique approach to 
study the raw computational effort needed to generate 
such entanglement. 



2 



From an experimental standpoint, a QCA has a signif- 
icant advantage over a QC because individual qubits in 
the lattice do not need to be separately addressed, since 
uniform rules are applied in parallel across the lattice. In 
such an implementation, applying uniform fields over the 
entire system helps to eliminate error resulting from cross 
talk on neighboring qubits due to imperfectly aligned 
control fields. Some specific physical systems have been 
proposed as candidates for QCA including quantum dot 
arrays [ll| and endohedral fullerenes 12]. Throughout 
the development of the general QCA formalism in this 
paper, we provide specific examples of possible experi- 
mental implementations in order to emphasize the rele- 
vance of the QCA approach to present day technologies. 

In Sec. II we introduce the formalism for QCA and 
show how to construct arbitrary three cell neighborhood 
rules using homogeneous pairwise interactions and sin- 
gle qubit gates. We show how to transport quantum 
information with QCA in Sec. Ill and demonstrate op- 
timal sequences to swap two distant quantum states and 
to prepare three types of entangled states. In Sec. IV 
we explore general properties of entanglement dynam- 
ics with QCA. The dynamics of multi-spin entanglement 
are measured by a function linearly related to the pu- 
rity of the single qubits averaged over the lattice. This 
measure has the advantage of being observable in a phys- 
ical system supporting the QCA architecture. In Sec. V 
we extend the theory to open systems and demonstrate 
how more general non-unitary rules can be implemented 
in the QCA paradigm using measurement and quantum 
feedback. It is shown that for a particular mixing of 
a non-unitary rule with a unitary rule, entanglement is 
generated across a spin chain where there is none for 
purely unitary evolution. This is an example of envi- 
ronment assisted entanglement generation. Finally, we 
present conclusions and open questions in Sec. VI. 



II. FORMALISM 

A. Simulating QCA rules 

Consider a ID array of n lattice sites occupied by 
qubits ordered to n — 1. We define a radius r QCA 
as one that changes the state of a qubit at site j de- 
pendent on the states of the qubits in the neighborhood 
[j ~ r -> j + r ] - Given a system with nearest neighbor in- 
teractions, the simplest unitary QCA rule has r = 1 
describing a unitary operator applied over a three cell 
neighborhood (j — 1, j, j + 1): 

MOoo, 1x01,^10,^11) = |00)(00| ®ix 00 + |01)(01| ®ix i 
+ |10)(10| (8)wio + |11>(11| <8)txii. 

(1) 

where \ab)(ab\ <g> means update the qubit at site j 
with the unitary u a b if the qubit at site j — 1 is in state 
I a) and the qubit at site j + 1 is in state \b). In classical 
CA the local update rule M can be applied in parallel 



to all cells. To do so requires that a separate register 
store the current state of the lattice so the previous state 
of the neighbors is known before the cells are updated 
in parallel. For instance, radius 1 CA rules could be 
implemented by copying the current state, updating the 
even ordered cells on the original and the odd ordered 
cells on the copy, and splicing the updated cells together. 
By the no cloning theorem 13], non-orthogonal quantum 
states cannot be copied so this is not possible for QCA. 
However, the update can be divided into two stages: first 
update all the even qubits with rule M, next update all 
odd qubits. This rule is denoted a Block partitioned QCA 
(BQCA) and guarantees that at each stage the operators 
commute and thus can be implemented in parallel Q. 

We show that any BQCA can be simulated with a lat- 
tice of even order constructed with an alternating array of 
two distinguishable species ABABAB . . . that are glob- 
ally addressable and interact via the Ising interaction. 
In deriving the construction of QCA rules we initially 
assume periodic boundary conditions (n + j = j). The 
simulation is shown to be easily adapted to a lattice with 
fixed boundaries. 

The general pairwise interaction Hamiltonian across a 
ID lattice is 

= E E ( 2 ) 

j=0 a ,/3=0 

where the Pauli operators are labelled {cr a } = 
{l,<j x ,<jy,<j z } and the g J a At) are, possibly time depen- 
dent, coupling strengths (H = 1). It is straightforward 
to show that in order for the Hamiltonian to commute 
over all nearest neighbor pairs with periodic boundary 
conditions, it is restricted to the form 

n— 1 n — 1 

Hi(t) = Es^K ® + E ( 3 ) 

3=0 3=0 

where ai = a • r« defines the local Bloch vector at site 
j. We identify the local basis of each qubit along this 
Bloch vector meaning a J z = cri . . The second sum in Eq. 
prefers to single qubit free Hamiltonians. Note that in 
order to satisfy the periodic boundary conditions, n must 
be even. 

The system dynamics can be controlled in a non-trivial 
way with limited addressability by assuming a ID lattice 
constructed two species of qubits A and B arranged in 
antiferromagnetic order. Here the species may have dis- 
tinguishable two level energy spacings, h? = h A ^ for 
j even(odd), meaning the species are addressable in fre- 
quency allowing even or odd ordered qubits to interact 
in parallel with an external field. The two species could 
also correspond to disjoint two dimensional subspaces of 
the same four dimension system. In either case, a general 
control Hamiltonian that performs single qubit rotations 
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on the two species is written: 

n/2-l n/2-1 

H c (t) = n A (t) • ° 2J + fi B (t) • * 2i+1 - ( 4 ) 

j=0 j=0 

The total Hamiltonian acting on the system is H = 
Hi(t) + Hc(t). For simplicity, we assume an isotropic 
pairwise interaction gi (t) = g(t) corresponding to the 
Ising interaction, and transform H to the appropriate 
rotating frame so that the total Hamiltonian becomes: 

n-l 

H\t) = H'jW+H'c® = g(t) oi®<d +1 +H'c(*)- (5) 

3=0 

Discrete time dynamics describing cellular automata can 
be implemented with continuous dynamics by first evolv- 
ing the system with the interaction Hamiltonian followed 
by evolution by the control that performs simultaneous 
single qubit gates on either or both species. This can 
be realized with a fixed Ising interaction punctuated by 
"hard" control pulses as is done in NMR pulse sequenc- 
ing , or the physical system may allow the pairwise 
couplings to be turned off during the single qubit gates. 
In any case, because the interaction Hamiltonian com- 
mutes with itself at all times, the unitary corresponding 
to coupled evolution can be written U(t) = e~ 1 ^ H i(t) dt 
and the single qubit gates are generated by He- 

The simplest non trivial BQCA rule obtainable in the 
two species architecture is described by the following gate 
sequence: 

M(l,u,u,u 2 ) = e- i ^ a ^U([y/(2g)])e^ az e i ^, (6) 

where u = e ll(Tft is an arbitrary element of 
SU(2) written as a rotation about the Bloch vec- 
tor ft = (sin # cos (/>, sin 6 1 cos 0, cos#), and m = 
(sin#/2cos(/>, sin#/2cos</>, cos 0/2). The evolution time 
in brackets, viz. U(t = [x]), means add the appropriate 
multiple of ir/\g\ to the quantity x to make it positive. 
Henceforth, we assume g > 0. All qubits of species A(B) 
are updated by the rule M A ^ when the single qubit 
gates act on that species. A single step of the BQCA is 
defined as the sequence: M = M A M B which updates all 
cells, and the BQCA rule iterated t times on the initial 
state |^(0)) generates the state \i/;(t)) = [M]*|^(0)). 

The particular rule Eq. ((SJ) is left /right symmetric 
with the interpretation that the center qubit is rotated 
by an amount proportional to the sum of the qubit 
values of the neighbors. Note that M A (u 2 , u, u, 1) = 
erf M A (1, u, u, u 2 )a B . Another elementary rule is: 

M(l,u,u,l) = -e- i t a -[e- i f^/7([7/(4^)])e i (f-2)^ 
e i2 i (J yU([-37r/(4g)})e i ^ (7z e- i2 i (Jz e i 2 ( 'y 
U{[-tt / {Ag)))e l ^ a y e^ az e^-i^ 
e i f^ ze i ^ a yU([y/(4g)])e- i i az ]e i 2^^ 

(7) 

where the unitary u and the Bloch vectors (n, m) are de- 
fined as above. Combining rules|H|and[3we can construct 



all symmetric QCA rules: 

M(^oo, u i = ^10,^11) M(v 2 , v, v, l)M(l,w,w, 1) 

M(l, u, u,u 2 ), 

(8) 

u 1/2 1/2 1 -1/2 -1/2 a 

where v = u 00 ,u = u xl , and w = u 0Q uoi^n • A 
maximum of six pairwise interactions U interspersed by 
single qubit gates is sufficient to simulate the symmetric 
rules. 

Asymmetric rules can be constructed if the Ising in- 
teraction is allowed to have different coupling strengths 
between left-center and center-right pairs. The appropri- 
ate Hamiltonian is: 

Hasym = ^\g\t)ay®al^+g\t)al^®al^). 

3=0 

(9) 

This asymmetry can be built into the system as is sug- 
gested, for instance, in proposals to implement quan- 
tum computation in optical lattices 15]. Here atoms 
are trapped in a 3D periodic potential created by stand- 
ing waves of interfering laser beams and prepared with 
one atom per potential well. An antiferromagnetic or- 
dering of atomic species can be created along one dimen- 
sion, and by appropriate tuning of the laser parameters, 
wells can be joined along this dimension such than each 
atoms interacts with its left or right neighbor. By choos- 
ing different interaction strengths (or times) between the 
center- left and right neighbors, H asyra can be simulated. 

Given the ability to implement H asym , the following 
rule can be generated: 

M(l,u-\u,l) = e-t^Uiit])?***, (10) 

where the time and couplings satisfy 

/ V(*)+ff 2 (*))d* = (mod27r), 1 } 

and (u, m, ft) are as above. A general QCA rule can be 
constructed from the above elementary rules: 

M(^oo,^oi,^io,^n) = M(l,x _1 ,x, l)M(l,x,x, 1) 
M(v 2 , v,v,l)M(l,w,w,l) 
M(l,u,u,u 2 ), 

(12) 

-1 1/2 / -1/2 -1/2 , 

where v = u 00 ,x = y uiou 01 ,w = u 00 uoi^n , and 

1/2 

u = . A maximum of 11 sequences of left/right shifts 
punctuated by single qubit gates are sufficient to imple- 
ment an arbitrary 3 cell BQCA rule M, although it is 
uncertain whether this is optimal. 

The present construction of M(izoo> ^01 > ^10? ^11) with 
Uij G SU(2) is only a subset of the most gen- 
eral rule having three additional relative phases, 
M(iioo, e^iioi, e^ 2 iXio, e^ 3 wn). One relative phase can 
be fixed by applying a z rotation, e^^, to the neigh- 
boring qubits after an update. Also, for each unitary 
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u = e* 7cr ™ the replacement 7 — > 7 + 7r introduces a 
sign change. A second relative phase is fixed by apply- 
ing different controlled phase gates between a qubit and 
each neighbor (as is possible with the interaction H asym ). 
Generating a third relative phase requires a direct inter- 
action between left and right neighbors, not possible in 
the ID architecture with only nearest neighbor connec- 
tivity and two species addressability. 

B. Boundary Conditions 

In the above treatment we have assumed periodic 
boundary conditions, in which case the the BQCA rules 
can be implemented uniformly with only global address- 
ability of two species. In practice, a system coupled by 
the Ising interaction with periodic boundaries could be 
realized by a custom designed ring molecule with alter- 
nating atomic species, or perhaps with trapped atoms in 
a ring type cavity. In most experimental situations it will 
be easier to construct a linear system with boundaries. 
Consider an open ID spin chain labeled from left to right 
by the integers to n — 1, where as above n is assumed 
even. One can simulate evolution where each cell is up- 
dated according to neighborhood values by introducing 
fictitious boundaries on the left and right ends with fixed 
values: G {0,1}. This is accomplished by append- 

ing appropriate single qubit gates to the ends of the chain 
after each instance of U (t) in the above rules. For a se- 
quence updating even ordered (A) species, append U (t) 
with e - *( -1 ) Lgta z, and for a sequence updating odd or- 
dered (B) species, append with e - ^ -1 ) R 9 ta z . I n this 
way, addressability at the ends (or at least the ability to 
introduce energy shifts at the ends) is sufficient to simu- 
late BQCA rules over ID systems with boundaries. 

C. Universality 

We have identified a finite set of rules to construct a 
class of radius 1, two state BQCA. We now discuss some 
issues regarding the universality of this class of quantum 
cellular automata. By universality we refer to the ability 
to emulate other computations, in particular other QCA 
and quantum computers, in an efficient manner. A dis- 
tinguishing feature of computation with classical cellular 
automata is that CA have minimal time complexity in 
that the same rule is applied to the data register at each 
iteration. This is in contrast to the conventional comput- 
ers that use a complex sequence of logic gates over the 
period of computation. It has been shown that a radius 
1, two state classical CA rule, designated rule 110, is uni- 
versal in the sense that by appropriate choice of initial 
state it can emulate any other CA as well as a turing 
machine Q. It should be emphasized that this rule up- 
dates all cells synchronously. It is not obvious that by 
appropriate choice of initial state, a single BQCA rule 
would be universal in the same way. It can be shown, 



however, that a sequence of rules can simulate a quan- 
tum computer with only linear cost in space and time re- 
sources. This may violate the spirit of using a single rule 
to generate complex dynamics but demonstrates that the 
underlying physical architecture for QCA supports uni- 
versal computation. The ability to simulate a QC follows 
from the work of Benjamin who has shown [^.Hq that an 
open ID lattice composed of an alternating array of two 
species of qubits can be used for quantum computation. 
The only architectural requirements are global address- 
ability of the species and addressable boundaries. He 
proposes two protocols with different physical assump- 
tions. 

The first protocol [3 assumes that unitaries SJ can be 
implemented in parallel, meaning "apply the unitary u 
to species S if the field value is equal to /". The field 
value is defined as the number of O's minus the num- 
ber of l's in the neighborhood of each qubit, therefore 
/ G {—2,0,2} inside the lattice and / G { — 1,1} at the 
boundaries. BQCA rules provide an explicit construc- 
tion of these unitaries, viz. the sequence S^SqS^ acting 
on interior spins is exactly simulated by the symmetric 
BQCA rule M s (u, v, v, w). In this proposal the boundary 
spins need to be addressable in order to load information 
into the lattice, but the only operations needed are the 
bit flips, S°J\ When the entire lattice is coupled via 
the Ising interaction, this is achievable by dynamically 
decoupling the boundary spins from the rest of the lat- 
tice using standard refocusing techniques 14]. Explicit 
pulse sequences to perform computation with endohedral 
fullerenes in a QCA architecture have been recently pro- 
posed in P3. 

In the second protocol 16] it is not necessary to apply 
unitaries that depend on the total field value, however, 
it is assumed that the Hamiltonian between left /center 
neighbors, H AB , can be turned "off' while the cen- 
ter/right Hamiltonian, H BA , is turned "on", and vice 
versa. This is akin to the physical requirement for asym- 
metric BQCA rules and may be more difficult to engineer 
in a given system. 

The ability to map BQCA to Benjamin's model of com- 
putation resolves a question about whether BQCA are 
universal with respect to the ability to efficiently simu- 
late other quantum cellular automata. Watrous p has 
shown that ID-partitioned QCA can be simulated by a 
quantum turing machine QTM) with only linear slow- 
down. ID-partitioned QCA are a restricted class of 1D- 
QCA in which each cell is partitioned into three subcells 
and the rule updates the cells by permuting subcells of 
neighboring cells and operating on the new cells in par- 
allel with quasi-local unitary operations, van Dam [l7| 
extended this result to prove that quantum gate cellu- 
lar automata (QGCA) can simulate any unitary QCA 
with only a polynomial slowdown. QGCA evolve by a 
repeated sequence of two steps: one step acts to per- 
mute the basis states within a certain neighborhood, and 
the second step applies parallel quasi-local gates over the 
neighborhood. Not all of the quasi-local gates in the 
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QGCA model can be implemented with pairwise interac- 
tions and two species addressability. However, one can 
use Benjamin's protocols to show that with a properly 
prepared initial state, a sequence of homogeneous update 
rules can simulate the QGCA model with only a linear 
cost in space and time resources. Because BQCA rules 
are sufficient to implement Benjamin's model, BQCA are 
also universal in this respect. 



III. INFORMATION TRANSPORT 

The discrete time process corresponding to QCA evolu- 
tion is a useful way to study information flow in quantum 
systems. For classical CA the maximum speed of infor- 
mation flow, c max , is 1 cell per update, which defines a 
light cone for information propagation. This can be real- 
ized, for instance, by beginning in the state . . . 010 ... 
and evolving with Wolfram's rule 254 4]. This rule up- 
dates the center cell in a three cell neighborhood, map- 
ping each cell to a 1 unless its left and right neighbors 
are in state 0,0. Evolving the initial state will cause the 
string of l's to grow by one cell on the left and right at 
each step. This does not fit into the QCA paradigm for 
two reasons. Firstly, the local rule is not unitary, e.g. 
both strings "110" and "111" are mapped to "111". Sec- 
ondly, the applied rule updates all cells simultaneously, 
not in a block partitioned manner. There are locally re- 
versible CA rules that spread information at speed c max 
(such as rule 150), but is there a unitary BQCA rule that 
can saturate the speed limit? The answer is affirmative 
as is shown below. 

We consider an n cell register initialized in the state 
| l)o ® |0)i... n -i with left and right boundary conditions 
(TtL — 0, Y,r = 0). The approach is to map the soli- 
tary 1 into a two cell unit which then propagates 2 
sites per update and is decoded into a single cell at the 
other boundary. The BQCA sequence to achieve this is: 
^^[M^e^^^-^^^-^Jf/ 2 . The total time 
to transport the information over n — 1 cells is (from Eq. 
EJ T n = nir/ (4<7) + (n+l)t s , where t s is the time to imple- 
ment a single qubit gate. If we assume that single qubit 
gates can be implemented on a time scale much faster 
than the many body interaction, t s <C 7r/(4g), then the 
information speed is c = (n — 1)/T ~ Ag{n — l)/(7rn). 
The connection is made to classical CA's by noting that 
nontrivial rules use conditional bit flips , which, accord- 
ing to Eq. El are implemented in a time t = 7r/(4g). 
Thus the maximum speed of information flow is c max = 
1/t = 4g/7r. The encoding and decoding consume a fixed 
amount of time but in the limit of large n, c max is ap- 
proached. 

The state could just as well have begun in a superposi- 
tion state \(j))o = a|0)o + /3| l)o ? m which case the BQCA 
sequence will transport the state to site n — 1 as shown 
in Fig. H It does so by first mapping the product state 
into a two particle entangled state, then shifting the ad- 
jacent l's two cells per update and finally mapping back 



to a product state at the other boundary. The entangle- 
ment present during transport is evident from the space 
time diagrams in Fig. ^ The first diagram shows the 
probability density for each spin to be in state |1) de- 
fined by Pi(pj) = Tr[\l)jj(l\pj], where p { = Tr not *(p) 
is the reduced state of the spin at site i of the global 
state p. The second diagram displays the reduced von 
Neumann entropy, defined by S(pi) = —Tr(pi\ogpi). 
Starting from a pure separable state and evolving uni- 
tarily, the residual mixedness of each spin in the ad- 
jacent pair results from mutual entanglement. We ex- 
plore in more detail the dynamics of entangled states 
below. By linearity, using the above BQCA sequence, 
any state po can be swapped through n — 1 lattice sites 
prepared in |0) in a time T n ~ mr/(4#). At the cost 
of one additional update on the B species, the states 
of two qubits on the ends of a chain can be swapped 
via the sequence: cr^cr^ 1 M B (1, e~ l % CTx , e~ t ^ CTx J e~ t7V(Tx ) 
[M(l,e- l f^,e-^^,e-^)] n / 2 , as shown in Fig. 

The QCA state transport time T n = n7r/(4g) is prov- 
ably 18] the minimal time to translate a quantum op- 
erator over n — 1 sites of a spin chain with pairwise 
interactions. This protocol should be compared to the 
total propagation time using "soliton operators" pro- 
posed in pjj. There a single qubit state is encoded 
into three qubits which then propagate through the chain 
and are decoded to one qubit at the end. Similar to the 
present proposal, the discrete evolution is generated via 
the Ising interaction punctuated by homogeneous single 
bit gates. The swapping time using the soliton opera- 
tors is T so i n = (n + l)7r/(4g), slightly longer than the 
present method. Both pulse sequences require address- 
ability at the boundaries but the simpler QCA sequence 
has the additional requirement of an architecture that 
supports an alternating array of two species. It should 
be emphasized that neither of these methods are true 
swap sequences, in the sense that any quantum informa- 
tion encoded in the intervening cells will be disturbed 
during the sequence. They may be useful in quantum ar- 
chitectures where quantum "memory" is stored in qubits 
spatially separated from each other by "bus" qubits ini- 
tialized to the state |0) that act as conduits for quantum 
information. Architectures with this kind of sparsely dis- 
tributed memory avoid correlated errors induced by the 
environment and can make the system amenable to quan- 
tum error correction. 

BQCA rules also can readily be constructed to 
distribute entanglement. Consider the creation of an 
entangled pair of qubits at the boundaries of an open 
chain of size n > 4. Choosing boundary conditions 
(Ex =0, E# = 0), we begin with a single qubit "seeded" 
to the superposition state l/\/2(|0) + |1)) near the 
middle of the chain with all other spins initialized 
to the state |0) and apply a QCA sequence to create 
the maximally entangled pair described by the state 
1/V2(|0) <8> |0) n _i + |1) ® |l)n-i). The particular 
BQCA sequence and optimal location of the seed spin 
will depend on the size of the lattice n. We choose 
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probability Pi(pj) reduced entropy S(pj) 
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space space 



FIG. 1: Transporting a quantum state over an n — 
14 cell ID lattice via BQCA evolution by the rule 
M(l,e~ i ^ C7x ,e~ i ^ ax ,e~ i7T(Tx ). Space time diagrams are 
shown with cell number on the horizontal axis and time flow- 
ing downward. On the left is a history of cell site probability 
to be in state |1) and on the right, the reduced von Neumann 
entropy of each cell (Black=l,White=0). 
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FIG. 3: Generating entangled states beginning with a 
"seeded" qubit in the superposition state l/y/2(\0) + |1}). (a) 
A maximally entangled pair at the boundaries: l/v / 2(|0)o 
|0) n -i + | l)o |l)n-i). (b) The n-spin Greenberger-Horne- 
Zeilinger (GHZ) state: 1/V2(\0 . . . 0) + |1...1». 
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FIG. 2: Swapping the states of cells at the boundaries of an 
n — 14 cell lattice through intermediary cells initialized to 
10} • 



a convention that this spin be of the A species and 
closest to the middle of the chain. For n = 4/c, k G N, 
the seed spin is located at site n/2 and the update 
sequence is: e _i ^°M B (l, , e -*f °* , e -™°* ) 

[M{1, e-^ Gx , e _< 5 Gx , e~ i7r(Tx 

M(l,e _i f (T X ^ e -ifa x ^ e -i^cT X y Similarly, for n = Ak + 2, 
the seeded spin is located at site n/2 — 1 and the update 
sequence is: e"*?^ [M(l, e~ 1 ^ Gx , e~ l ^ (Jx , e ~ l7rcTx )] k 
M(l, e~ z % (Tx , e~ z % (Tx , e~ l ^ (Jx ). The sequence works by 
updating the state |0 . . . 010 . . .0) once with the rule 
M(l, e~ l2 5 ax , e~ l2 5 ax , e~ % ^ Gx \ creating a separated 
pair of adjacent spins in the state |1). These pairs 
then propagate outward under the same rule used to 
transport quantum information to the boundaries. An 
example for n = 14 is shown in Fig. EJa). The total 
time to produce an entangled pair at the boundaries 
is calculated, using Eqs. 16171 and again assuming t s is 
negligible, to be T n = (4 + n/2)7r/(4g). By a similar 
argument to the optimality of quantum state trans- 
port sequence, the present sequence for distributing 
entanglement is optimal within the QCA framework. 
However, if one is allowed to perform measurements 
as well as unitary evolution, then an entangled pair 
can be produced at the ends using "entanglement 
swapping". Given the same architecture under present 



consideration, it has been shown 3] that a maximally 
entangled pair can be swapped to the ends of a spin 
chain in a time, T = n/(2g) independent of the length. 
The entangled pair could then be used as a resource to 
teleport a quantum state from one end of the chain to 
the other. Naturally, any protocol to deterministically 
distribute quantum states must preserve causality and 
is fundamentally limited by the speed of light which 
enters into the protocol through classical processing of 
measurement results over the length of the chain. 

Multi-particle entanglement can be constructed using a 
slight variation of the sequence for distributing entangled 
pairs. As above we assume boundary conditions (El = 
0, En = 0) with a seeded qubit in the superposition 
l/V2(\0) + |1)). An n-spin GHZ state 1/V2(|0 . . . 0) + 
| X ... X ) ) can be generated as follows. For n = 4/c, k G N, 
the seed spin is located at site n/2 and the update se- 
quence is: e"^" 1 )*^ [M(l,e- i ^ ax ,e- i ^ <TiB ,e~ i ' ir<Ta! )] k . 
Similarly, for n = Ak + 2, the "seeded" spin 
is located at site n/2 — 1 and the update se- 
quence is: e t(_1) H^M B (l, e _i t ° x , e~*t Gx , e - i7TCTx ) 
[M(l,e~ i 2 <T ',e~ i ^ <TiB ,e~ i ' ir<Tx )] k . An example for n = 14 
is shown in Fig. 0{b). The total time for the BQCA 
sequence is: T n = n7r/(Sg). 

Choosing other, more complicated, initial states can 
allow BQCA generation of many classes of multi-particle 
entanglement. For example, if we fix boundaries at 
(El = 0, E# = 0) and initialize the n cell lattice to 
= e~ % ^ ^ j=o a y |0 . . . 0), then the resultant state after 
one update by the rule: M(l, e~ % ^ Gz , e~ % ^ Gz , e~ l ^ az ) is 
characterized by many-particle quantum correlations. In 
fact it is equivalent, up to local unitaries, to the so called 
cluster state [20j : 

1 n — 1 

|*n clue) = (g)(|0> o <X« +1 + |1) ), (13) 

a=0 

with the convention <j n = 1. These states obtain max- 
imal reduced entropy of every spin. Most notably, they 
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have the property of maintaining persistency of entangle- 
ment between the remaining set of qubits when some are 
lost (depolarized, measured, etc.). They have exponen- 
tially large Schmidt number, namely any expansion of 
the state in terms of a product basis will require at least 
2 n / 2 terms. This is to be contrasted with the n-spin GHZ 
state which has Schmidt number 2. 



IV. ENTANGLEMENT DYNAMICS 

A. Quantifying Multi-Spin Entanglement 

Generally, QCA evolution can take a configuration of 
spins prepared in a product state to a number of different 
entangled states. In order to characterize the dynamics 
of entanglement, it would be helpful to have a single pa- 
rameter that quantifies the amount of multi-particle en- 
tanglement contained in a state at any given time step. 
A good measure of entanglement should capture the non- 
local nature of the quantum correlations of the spins and 
therefore should be a function on the state that is non in- 
creasing, on average, under local operations and classical 
communication. Because entanglement can be shared in 
different ways by different subsets (parties) of the spins 
in the lattice, this is no single functions that describes 
multi-partite entanglement. For the purposes of this pa- 
per we quantify the amount of multi-spin entanglement 
with a function on pure states of n qubits introduced in 
[2lj and expressible as: 

«» = 2^1 -1/n J2 Tr^. (14) 

The measure R is linearly related to the purity of the 
single qubits averaged over the lattice and satisfies two 
important properties. First, < R(\ip)) < 1, where 
R(\ip)) = iff is a product state, and R(\ip)) = 1 
for some entangled states. Second, R(ijj)) is invariant 
under local unitaries Uj. 

A significant advantage of this function over other pos- 
sible measures is that it can be observed in a straight- 
forward manner by measurement. This can be done by 
introducing a second, identical ID lattice and interact- 
ing the two lattices, bitwise, with a third addressable ID 
lattice that can be prepared and measured. The mea- 
surement requires only that each lattice be addressable 
but does not require addressability of cells within the 
lattice and is described in [22j. A deficiency of R as an 
entanglement measure is that it cannot distinguish sub- 
global entanglement. For example, in a n = 4 lattice, 
the product state of two maximally entangled Bell states 
and the 4-spin GHZ state both have R values equal to 
1. This should be kept in mind when quantifying the 
entangling capacity of BQCA rules as is done below. In 
principle, there are other measurements that can be car- 
ried out over a many spin system to distinguish one type 
of shared entanglement from another. 



B. Generating multi-spin entanglement 

In Sec. Ill, we considered some examples of BQCA 
rules that generate and distribute entangled states. In 
this section, instead of searching for BQCA sequences 
that generate particular entangled states, we explore 
some basic properties of the rules themselves by way of 
two examples. While these examples are not intended 
to simulate any particular physical system, they do illus- 
trate some universal behaviors of BQCA and indicate the 
computational power of simple rules applied over local, 
in this case 3 cell, neighborhoods. 

It would be beneficial if some predictive statements 
could be made about the behavior of QCA. We know 
that classical cellular automata have the property that a 
globally reversible rules follow closed evolution 4] . That 
is, any initial configuration will evolve back to itself after 
a characteristic period that depends on the rule and the 
configuration. The maximum period of evolution is the 
size of the configuration space, which for an n cell lattice 
with 2 states per cell is 2 n . Linear CA rules are those that 
satisfy the property that for an initial configuration that 
is a mixture of two configurations, u = av + bw : where 
a, b G R, the rule acts linearly on the inputs: Mu = 
aMv + bMw. If the periods of the configurations v and 
w under rule M are Tm(v) and Tm(w) respectively, then 
the period of u is Tm(u) = Icui(Tm (v) , Tm (w) ) • For an 
arbritary mixture of m configurations, u = Y^k=i a kVk, 
the period is Tm(u) = lcm({TM(vk})- 

We consider the BQCA rule Ml = 
M(l, e~ l ^ ax , e~ l % (Tx , e~ l % (Tx ) acting on an initial 
product state of spins with boundaries fixed at 
(Dl = 0, Y>r = 0). For the computational basis states 
with |0) on the interior qubits and |0) or |1) at the 
ends, the evolution is reversible as shown in Fig. 
These four initial states have characteristic periods 
1,11, and 13 and at no time is entanglement gener- 
ated by the rule. If the initial state is chosen as an 
evenly weighted superposition of these four basis states: 
|^(0)) = e~ i 5 (<T S+°"?" 1 )|0. . .0), then entanglement is 
generated by the rule because non-separable phases 
accumulate on the co-evolving basis states. The space 
time diagrams of probability density and reduced 
entropy are shown in Fig. EJl The evolution is periodic 
with a period given by T = Zcra(l,ll,13) = 143. The 
multi-qubit entanglement during the evolution is plotted 
in Fig. The entanglement never attains values above 
R(\ip(t))) = 0.6, and this is evident in the space time 
plot of reduced entropy which shows that at any given 
time step no more that 8 out of 10 spins are entangled 
with each other. 

The entanglement dynamics are dramatically different 
for the same system and initial state evolving under the 
BQCA rule M2 M(l, e _i * a * , °% e - ^). This 
rule rotates each spin by half the amount of the rule 
Ml, however, the quantum dynamics does not follow a 
simple composition rule, i.e. Ml ^ (M2) 2 . This is be- 
cause the underlying Hamiltonians that update the states 
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(a) P 1 ( Pj ) (b) P 1 ( Pj ) (c) P 1 ( Pj ) 




FIG. 4: Evolution of three computational basis states by the 
rule M(l,e~ i % (Tx ,e~^ CTa; ,e~ i ^ (Tx ) over an n = 10 lattice with 
boundaries fixed at |0). Shown are the probability density 
space time diagrams over one period for the initial states: 
(a) 1 10... 00), (b) |0...01), and (c) |10...01). The evolution 
of the state |0 . . . 0) is trivial. At no time is entanglement 
generated. 



of species A and B are non-commuting so that correla- 
tions build up at rates that are not linearly related. The 
space time diagrams of probability density and reduced 
entropy plotted in Fig. show that after three steps, 
correlations spread throughout the lattice. No periodic- 
ity is evident, and after roughly 20 steps, the multi-qubit 
entanglement, plotted in Fig. saturates at a value of 
approximately R(\i/j(t))) = 0.9 with small fluctuations. 
One way to discern whether the rule is generating many 
different classes of entangled states during the evolution 
is to examine temporal variation of Schmidt numbers of 
the state \ijj(t)) over the set of all 2 n_1 bipartite divisions 
of the n lattice qubits. The Schmidt numbers are invari- 
ant under local unitary operations and under a bipartite 
division of k and n — k qubits, their range is the integers 
in the interval [0, mm{2 fc , 2 n ~ k }]. We have calculated 
the history of Schmidt numbers over the evolution pe- 
riod and find that the rule generates a large number of 
different classes of entangled states. This demonstrates 
that the rule M2 explores a larger volume of the Hilbert 
space of pure states in Hf n than does rule Ml for the 
given initial state. 

It is an open question under what BQCA rules and ini- 
tial states is the set of states generated during evolution 
dense on the Hilbert space of pure states in TL® 71 . One 
might expect that rules that rotate the updated spins by 
an angle that is an irrational multiple of tt would accom- 
plish this for a large class of initial states. 

V. NON-UNITARY RULES 

A. Formulation 

Up to this point we have described how to implement 
a class of unitary BQCA rules. In general one would like 
to have a prescription for implementing non-unitary rules 



Pi(Pj) S( Pj ) PM) S( Pj ) 




space space space space 



FIG. 5: Entanglement dynamics visualized by the space time 
histories of the evolution of a chain of 10 spins by two BQCA 
rules. The boundaries are fixed at |0) and the initial state is 
the same for both rules with all qubits initialized to 0) except 
for the qubits at sites and n — 1 each in the superposition 
state 1/V2(|0) + |1». (a) Rule M(l, e -4 *^, eT^** , e~ i ^ C7x ). 
(b) Rule M(l,e-^ ax ,e-^ ax ,e~^ ax ). 
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FIG. 6: Multi-qubit entanglement generated during the 
BQCA evolution plotted in Fig.|5] Entanglement of the global 
state is plotted for the rules M(l, e~ i ^ (7x , e~ i ^ (Tx , e~ i ^ (7x ) 
(triangles) and M(l, e ~ i ^ (Tx , e~ i ^ (7x , e~ i ^ (Tx ) (boxes). 



as well. Of the 256 Wolfram rules for radius 1, two state 
classical CA, only sixteen are locally invertible. They are 
given by: 

U jtkt i tm = |00)(00|® + |01)(01|<8>(<7 a! ) fc ( . 

+ |io)(io|0Ky + |ii)(ii|0(^r, ^ } 

where {j, fc, Z,m} G {0,1}. Among classical CA, all the 
unitary rules generate rather simple behavior compared 
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to the complex dynamics generated from some of the 
other, non-unitary, rules 4]. For example, rule 110 is 
described by the following update table: 

oho. in HQ 101 100 on 010 001 000 
110 1110 

which is not one-to-one because both states 011 and 001 
map to 011. As mentioned above, rule 110 is univer- 
sal and we might expect that for quantum cellular au- 
tomata there are interesting dynamics to be explored 
when the register is no longer a closed system obeying 
unitary dynamics but interacts with an environment in 
an irreversible way. 

Non-unitary rules correspond to completely positive 
maps applied to a system dependent on the state of the 
neighborhood. A general, completely positive map on a 
quantum state p can be written as a superoperator in 
the Krauss representation 23] as: S(p) = Y^=i F^pF^ 
where the total number of effects is k and trace 
preservation of the state is ensured by the constraint 
SjUl F^Ffi = 1. In the QCA context, the effects acting 
on a three cell neighborhood (j — 1, j, j + 1) are a sum 
of actions on qubit j induced by orthogonal states of the 
qubits at sites j — I and j + 1. The superoperator can be 
written as the composition 

S J (p) = S™oS° 1 oS}°oS}\p), (16) 

where 

k a b 

Sf(p) = \ab)(ab\ ®^2f« b pff \ab)(ab\. (17) 

Here, fc a & denotes the number of effects that act on the 
updated qubit j when the neighborhood is in the state 
\ab)j-ij+i. The single qubit superoperators are trace 

preserving, i.e. J2^=i f^f^ = 1- As w ^ n the unitary 
maps, the maps Sj and Sj+2 commute, so qubits at every 
other site can be updated in parallel. We denote a total 
BQCA update sequence from time t to t+ 1 by: p(t+ 1) = 
%{p{i)) = % A o% B (p(t)), where 

§ A (p) = SooS 2 o---oS n - 2 (p), n «x 
$ B (p) = SioS 3 o...oS n _i(p). [ ^ 

As an example, the CA rule 110 updating the state of 
a qubit at site j is written 

R) w {p) = F 1 {j)pF\{j) + F 2 (j)pF}(j), (19) 

where 

F{ = |00)(00|<g)F' + |10)(10|(g)F 

+|11)(11|®^ + |01)(01|®|1)^(1|, (20) 
Fi = |01)(01|® 11)^(01. 

The rule can be decomposed into unitary and non- unitary 
BQCA rules as: 

Rf°{p) = Sf (M(l, 1, 1, 4)pM(l, 1, 1, 4)), (21) 



where 

S^( P ) = |01)(01|® (11)^(1^11)^(11 

+ |i)^(oho)^(i|)®|oi)(oi|). [22) 

When the neighborhood is in state |01), rule 110 has the 
effect of an amplitude damping channel on qubit j, i.e. 
it maps a mixed state to a pure state |1). 

B. Simulation 

In this section we demonstrate how to implement non- 
unitary rules within a QCA architecture. A non-unitary 
map on a quantum system residing in a Hilbert space TL S 
can be thought of as open system dynamics that arise 
from unitary operation in the combined space Ti s <S> Ti e 
of the system and some environment, followed by tracing 
over the environmental degrees of freedom. Any super- 
operator on a system of dimension d can be realized in 
this fashion with an environment of dimension at most 
d 2 ; meaning that the maximal number of effects in a su- 
peroperator expansion is k = d 2 [24| . 

Generally, implementing control over a combined sys- 
tem and environment of this size is difficult, however, it 
has been shown that by using measurement and feed- 
back, a single qubit environment is sufficient to simu- 
late open system dynamics |25| . The simulation works 
by coupling a single qubit, e, prepared in the state 
|+) e = 1/V2(|0) e + |l) e ) to the system, 8, via a Hamil- 
tonian H sym = 7P cr z , where P is a projector onto a 
pure state in s. The corresponding unitary operation is 
U(t) = e~ iHs y mt = cos(jtP) 1- ismfrtP) a z . By 
suitable averaging techniques, namely conjugating short 
time evolution, [/(At), where At <C t : with unitary oper- 
ations on 5, the projector P can be transformed into any 
positive, unit trace, operator P. Finally, e is measured in 
the a x basis and the result is fedback to s with one of two 
unitaries Uo,Ui, depending on the measurement result. 
This process will implement any superoperator described 
by two effects: Fo = Uocos^tP) and F\ = Uism(jtP). 
Maps with more than two effects can be simulated by 
repeated cycles of measurement and feedback. 

In the QCA context, we want to activate non-unitary 
dynamics on a qubit s dependent on the neighbor- 
hood state. This can be accomplished by first "turn- 
ing on" a control qubit c dependent on the neighbor- 
hood state, then implementing the unitary, U c - syrn = 
e -*|i) cc (i|(g) J H r sy?T ,^ during the simulation steps described 
above. For instance, the superoperator Sg 1 (p) can be 
implemented by first entangling the register and the con- 
trol qubit (initialized to |0) c ) with the unitary BQCA rule 
M c (l, 1, e~ l ^ ax ^ 1). We choose the projector in H sym to 
be P = |l) ss (l| and simulate open system dynamics on 
s by evolving the combined system (c, e, s) with the uni- 
tary U c - syrn = e-^^^e^^^e^^^. The 
gate U c - S y m can be efficiently implemented using pair- 
wise interactions between system pairs (c, e) and (e, s) 
[22| . At the measurement stage, the classical result of 
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FIG. 7: Sequence of steps to implement non-unitary BQCA 
rules on a ID lattice, s, using qubits in a control lattice c and 
an "environment" lattice e. The lattice cells are assumed ad- 
dressable along the vertical but not the horizontal direction 
with grey(white) cells corresponding to qubit species A(B). 
Shown is a sequence sketched over a three cell section of the 
lattices for implementing a non-unitary rule on B species 
qubits dependent on the A species neighbors, (i). A product 
state in the computational basis of s is shown with the lattices 
c and s initialized to |0). (ii). The state of the A species qubits 
is fanned out to corresponding sites in e and c using CNOT 
gates, (hi). The unitary BQCA rule M B (uoo, uoi, uio, un) 
acts on c to activate the controls dependent on the neighbor- 
hood. Non- unitary evolution on s is simulated using inter- 
actions between the lattices, measurement on e and feedback 
on s. Afterward, the inverse of steps (ii-iii) disentangles the 
three lattices. 



the measurement on e should be ignored, equivalent to 
tracing over the environmental degrees of freedom. This 
can be accomidated by an environmental qubit that sinks 
information of the measurement result to a large reser- 
voir (such as an atom that emits spontaneous radia- 
tion). The conditional feedback can be implemented 
without knowledge of the measurement result using a 
controlled- unitary operation between e and s, namely: 
U fh = |0) ee (0| <g> U$ + |l) ee (l| U[. Finally, the control 
qubit c needs to be disentangled from the register with 
the rule M c (l, 1, e l 2°"*, 1), and the "environment" qubit 
reset to |+) e . 

Note that this protocol has the unwanted effect of ap- 
plying a unitary Uq to the qubit s regardless of the state 
of the neighborhood, because the feedback is only condi- 
tioned on the environmental state which is initialized to 
|0). This can be obviated by instead using the feedback 
gate: U' fh = |0) ee (0| 0l s + |l) ee (l| ®U$ a Uf , and after the 
control has been disentangled from s, applying a unitary 
BQCA on s. For example, after the simulation of S® 1 
with the gate U' fh , the unitary BQCA rule M s (l, J7 , 1, 1) 
will apply the necessary feedback. By the superoperator 
decomposition, EqITtjI any non-unitary BQCA rule can 
be implemented by at most four instances of the above 
open systems simulation using the appropriate unitary 
rule M c for each neighborhood value \ab). 

A possible implementation of non-unitary BQCA rules 
is shown in Fig. [7|for an architecture with three stacked 
ID lattices. The protocol for interacting registers in an 
optical lattice is outlined in [22j . 



(a) (b) 
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FIG. 8: System evolution by a mixture of unitary and non- 
unitary rules: S 01 (p,p). Results are shown for p = 1/2 
(boxes) and p = 1 (triangles) with connecting lines to guide 
the eye. (a) Mixedness of the system p(t). (b) Entangle- 
ment, quantified as the average tangle r over all spin pairs. 
Increased coupling to the environment increases the entangle- 
ment. 



C. Results 

We investigate the effect of adding decoherence to a 
unitary BQCA by mixing the rules R 110 and R 108 = 
M(l, 1, 1, crj). This is described by a one parameter map 
on the neighborhood (j — 1, j, j + 1) written as: 

Sj(p,P) = Sf(M(l, 1, 1, 4)PM(1, 1, 1, (23) 
where 

S?(P,P) = |01)(01|® ((11)^(11 

+ v / r Tr p|o> JJ -<o|)p(|i) JJ <i| + vT^io^od 

+ ^11)^(0^10)^(11) ®|01><01|). 

(24) 

In the case p — > 0(1), the rule approaches R 108 (R 110 ). 
We study the evolution of entanglement under this map 
when the initial state is chosen to be the superposition of 

all computational basis states: 1^(0)) = e"* 1 ^^ ^(O). 
Note that R 108 permutes computational basis states so 
that |^(0)) is an eigenstate of this rule and no entangle- 
ment between spins is generated 1 . However, when the 
system is subjected to non-unitary rules, entanglement 
can develop. This behavior is illustrated in Fig. [51 for an 
n = 6 lattice with boundaries fixed at (E^ = 0, = 0). 
Two global quantities of the spin chain are plotted; the 
mixedness of the state, 1 — Tr[p 2 ], and the entanglement. 
The entanglement over the multi-partite mixed state p 
is calculated by averaging the pairwise tangle over all 
spin pairs (i, j). The tangle [26| is a monotonic function 
on pure or mixed states of two qubits assuming the value 
for separable states and 1 for maximally entangled Bell 
states. It is defined as a function of the reduced state 
Pij of qubits i and j: = [max{\\ — A2 — A3 — A4}] 2 
where the A& are the square roots of the eigenvalues, in 



1 Because a x £ SU(2), rule 108 is not strictly within the class 
of implement able BQCA rules in ID. However one can use the 
rule M(l,l,l,e*"2°" a: ) instead and correct for the phase using a 
controlled phase gate in the non-unitary implementation stage. 
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decreasing order, of Pijpij. Here, p^j is the spin-flipped 
version of p^: pij = a y o- y p\ 2 o~ y ® (J y 

As the amount of coupling to the environment, quan- 
tified by p, increases, the amount of entanglement in the 
spin chain increases. Additionally, beyond a certain time, 
the mixedness decreases with increasing coupling. This 
latter fact is because rule 110 acts as a neighborhood de- 
pendent amplitude damping channel that decreases the 
reduced entropy of each spin state. In fact, for p = 1 the 
state of the system relaxes after a time t = 4 to dynam- 
ics with period 6 and constant mixedness. For the mixed 
rule case, p = 1/2, in contrast, the system is not driven 
to periodic evolution but on long time scales experiences 
small fluctuation in mixedness and entanglement. 

The development of entanglement during open system 
dynamics, absent in the closed system dynamics can be 
attributed to the fact that the non-unitary rule acts to 
dampen the state of a qubit only if the neighborhood is 
in the state |01). The implementation of this rule re- 
quires a three body interaction followed by single spin 
decoherence. The net effect is to project the state to one 
that has some fraction of entanglement. This is an ex- 
ample of environment assisted entanglement preparation. 
Another, well known, example of such a phenomenon is 
the relaxation of two independent radiating dipoles into 
a maximally entangled sub-radiant state |2jj when the 
dipoles are close enough together to see the same electro- 
magnetic field. The results here show that entanglement 
can develop even when the environment acts only on one 
member of a neighborhood of spins. 

It is possible that this effect could be measured in 
the laboratory as a signature of neighbor dependent en- 
vironmental coupling. For instance, consider an array 
of three atoms of two species A and B trapped inside 
a high Q cavity with order ABA and aligned perpen- 
dicular to the cavity axis. The species are assumed to 
have a distinguishable set of two ground state manifolds 
\0)a,b and \1)a,b with different resonant excitation fre- 
quencies hu A ( B ) = E e A(B) - ^oa(b) and hu 1A ( B ) = 
EeA(B) ~ EiA(B) and possibly different decay rates Ja,b- 
If a laser field at frequency uol illuminates all three atoms, 
dipoles will be excited with dominant dipole-dipole in- 
teractions Vdd acting pairwise. Assume that the uol is 
extremely far off resonant to the \0)a,b —> \ 6 )a,b transi- 
tion so that dipoles are excited only when atoms of both 
species are in state |1). The interactions will shift the 
energy levels of the two species so that the effective de- 
tunings of the field will be ~ uj l — uo\ b — 2V dd /h, 
and A A ~ uj l - uj x A - V dd /h. For V dd large enough 
and appropriate choice of laser frequency, |A#| < jb 
while |A^| > 7a so that the field is in resonance with 
the excited state of atom B but not for the A species. If 
the resonant cavity frequency is close to c^ob, then the 
B species atom will preferentially decay to state |0)#. 
This type of decay corresponds to the non-unitary rule 
S = erf S^Of P°x iP) a x , in analogy to Eq. El The 
strength of V dd will determine the amount of coupling p. 



VI. CONCLUSIONS 

In this paper, we have shown how to construct a uni- 
versal class of radius 1, two-state QCA that are block 
partitioned and particularly suited for implementation in 
systems with a naturally endowed lattice type structure. 
In Section II A we introduced a universal set of BQCA 
rules in Eq. 1, which is the formal analog of the subset 
of 16 unitary rules of Wolfram's 256 radius one rules for 
classical CA. We demonstrated how these general rules 
can be simulated in a spin lattice with Ising interactions 
in conjunction with single qubit rotations (applied in par- 
allel across the lattice). This is an important result show- 
ing the physical relevance of the BQCA to experimental 
systems. 

Another important result is that we have suggested a 
new approach to quantum information: use BQCA to 
explore the raw computational properties of a physical 
system, such as the transport of quantum information 
and the generation of long-range quantum correlations 
throughout the system. This approach should be viewed 
as complementary to the standard treatment of quantum 
information processing centered around the quantum cir- 
cuit model of QC. In Section III we presented several 
specific examples of how BQCA rules can be chosen to 
distribute specialized entangled states across the lattice, 
for example, to create a n-spin GHZ state, shown in Fig- 
ure 3b. In order to visualize the development and spread 
of entanglement, we have introduced the idea of plotting 
a space- time diagram of the reduced entropy at each cell. 
We then explored more general entanglement dynamics 
in Section IV, with a focus on finding rules that are ef- 
fective at generating multi-qubit entanglement. For this 
task, we utilized the multi-qubit entanglement measure 
R(\ip)) related to the average purity of the constituent 
qubit s. 

Perhaps our most innovative contribution to the study 
of QCA is that we have developed the formalism in Sec- 
tion V that extends BQCA to open quantum systems, 
which evolve according to non-unitary rules. The BQCA 
rule is represented by an appropriate set of Krauss op- 
erators acting on the system density matrix. This more 
generalized treatment can be thought of as including the 
effect of correlated noise in the quantum evolution. From 
a practical standpoint, this extension is important for ex- 
ploring the effect of a broader class of errors than is typi- 
cally treated in the theory of error correction. On a more 
fundamental level, the non-unitary BQCA can be used as 
a testbed for exploring the interplay between quantum 
and classical computation-between quantum and classi- 
cal correlations in a discrete dynamical system. As a 
concrete example exploring this notion, in Section V we 
presented simulation results of a non-unitary BQCA rule 
that can be tuned continuously from a purely open irre- 
versible evolution using rule 110 (implemented in a block 
partitioned manner) to a purely closed unitary evolution 
using the quantum analog of rule 108. Remarkably, we 
uncovered the intriguing result that for a particular ini- 
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tial state entanglement generation is optimized by tuning 
the mixing parameter to a finite value, which suggests the 
possibility of environment assisted entanglement genera- 
tion. 

The are several outstanding issues regarding QCA that 
warrant future research. Can QCA be used to simulate 
complex classical or quantum dynamics? The study of 
QLG algorithms demonstrates that there are non-trivial 
classical dynamics of a single particle can be studied us- 
ing a lattice of spins. In these algorithms, the entan- 
glement generated is typically limited to local neighbor- 
hoods. It is worth while investigating whether large scale 
global entanglement generated by QCA can be used to 
advantage, perhaps to study properties of classical CA 
with quantum parallelism. Also the preliminary study 
of non-unitary QCA here suggests that in some cases 
open systems dynamics may be a more efficient way to 
navigate through Hilbert space than purely unitary dy- 
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